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Abstract 

We investigate the thermalization process of the Universe after inflation to determine 
the evolution of the effective temperature. The time scale of thermalization is found to 
be so long that it delays the evolution of the effective temperature, and the resulting 
maximal temperature of the Universe can be significantly low^er than the one obtained 
in the literature. Our results clarify the finite density corrections to the effective poten¬ 
tial of a scalar field and also processes of heavy particle production. In particular, we 
find that the maximum temperature of the Universe may be at most electroweak scale 
if the reheating temperature is as low as 6{T) MeV, which implies that the electroweak 
symmetry may be marginally restored. In addition, it is noticeable that the dark mat¬ 
ter may not be produced from thermal plasma in such a low reheating scenario, since 
the maximum temperature can be smaller than the conventional estimation by five or¬ 
ders of magnitude. We also give implications to the Peccei-Quinn mechanism and the 
Affleck-Dine baryogenesis. 
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Introduction and Summary 


1.1 Introduction 

Slow roll inflationary scenarios are successful in light of the solution to the horizon problem, 
flatness problem, and the origin of the large scale structure. Inflation is usually driven by 
a finite energy density of a slowly rolling scalar field, called inflaton. After the slow roll 
conditions fail, inflaton starts to oscillate around its potential minimum and the energy 
density of the Universe is then dominated by that of the oscillating inflaton. In order to 
proceed the big bang nucleosynthesis, the inflaton has to decay into radiation, which is 
referred to as reheating [ i, 2] . 

While the inflaton elegantly solves the above problems and provides the primordial den¬ 
sity perturbations, there still remains unsolved cosmological issues; dark matter and baryon 
asymmetry of the Universe. Scalar fields other than the inflaton may play essential roles in 
solving these remaining problems. For instance, dark matter can be explained by a pseudo- 
NG boson, called axion [3], which is associated with the spontaneously symmetry breaking 
(SSB) of PQ symmetry triggered by a PQ-charged scalar held [4, 5] . In supersymmetric 
(SUSY) theories, the Affleck-Dine mechanism can explain the origin of baryon asymmetry 
by using a baryonic scalar held [6, 7]. Generally scalar fields should have interaction terms 
with radiation in order to successfully lose their energy, otherwise they tend to dominate 
the energy density of the Universe. This very interaction makes the dynamics of scalar fields 
non-trivial. 

One of prominent effects caused by the interaction between the scalar fields and ra¬ 
diation is modification of their effective potential.*^ As the thermal plasma grows after 
inflation, it drastically affects the dynamics of scalar fields. At a sufficiently high temper¬ 
ature, for instance, thermal effects may induce a positive thermal mass for a scalar held 
and make it stay at the point of the potential where particles coupled to the scalar held 
remain massless. When the scalar held is responsible to the SSB of some symmetry, like the 
Higgs boson or a PQ breaking scalar held, the symmetry is restored at that high tempera¬ 
ture [14, 15, 16, 17, 18]. Then, a phase transition occurs when the thermal mass decreases 
down to the zero-temperature mass of the boson held. At the phase transition, topological 
defects, such as cosmic strings and domain walls, may form and affect the evolution of 
the Universe. In QCD axion models, there may be the domain wall problem when the PQ 
symmetry is broken after inflation [19, 20], while there is the isocurvature problem when 
it is broken before inflation [21, 22, 23]. The abundance of axion DM also depends on 
whether the PQ symmetry is broken before or after inflation. In SUSY theories, the Affleck- 
Dine mechanism can generate baryon asymmetry by using a scalar held carrying a nonzero 
baryon charge [6, 7]. The amount of baryon asymmetry is sometimes affected by the finite 
temperature effects [7, 24, 25, 26, 27]. These examples show that the thermal effects on 
boson fields drastically affect the evolutions of their dynamics and change their predictions. 
Therefore, we should clarify the evolution of the temperature after inflation. 

One might assume that the reheating occurs instantaneously and regard the reheating 
temperature as the maximum temperature of the Universe. More carefully one may solve 
the Boltzmann equation for the inflaton and radiation to include the production of radiation 
before the complete decay of inflaton [28, 29]. Even in this method, “instantaneous ther- 

For other relevant effects of background thermal plasma on scalar condensates in the Universe, see 
Refs. [8, 9, 10, 11, 12, 13] for instance. 
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malization” of radiation is implicitly assumed. In Ref. [30], however, we have pointed out 
that the time scale of thermalization of radiation is finite and affects the reheating process 
of the Universe after inflation.*^ 

In this paper, we further investigate the thermalization process of inflaton decay prod¬ 
ucts and provides the evolution of an effective temperature which describes the strength of 
finite temperature effects. We find that the maximal temperature of the Universe after infla¬ 
tion is much smaller than that obtained by the “instantaneous thermalization” assumption, 
by extending the analysis given in Ref. [30] to the era before the soft sector is thermalized. 
The discrepancy becomes larger for smaller decay rate of inflaton and can be as large as 
about five order of magnitude. 

Our results clarify the finite density corrections to the effective potential of the scalar 
field such as so called the thermal mass and the thermal log potential; which are essential 
so as to determine the beginning time of oscillation of the scalar field and also the condi¬ 
tion of the symmetry restoration by finite density corrections. For example, we find that 
the maximum temperature of the Universe may be at most ~ lOOGeV if the reheating tem¬ 
perature is as low as ^(1) MeV. In this case, the electroweak phase transition is marginally 
restored after inflation. We also find that the reheating temperature may have to be two 
order of magnitude lower than the PQ breaking scale in order to avoid the restoration of PQ 
symmetry. In addition, we investigate the effect on the AD field, which generates baryon 
asymmetry by the Affleck-Dine mechanism. We find that the result is consistent with the 
case that one assumes “instantaneous thermalization” when the VEV of the AD field is larger 
than the mass of inflaton. 

In order to study the finite density corrections to the effective potential, we have to 
discuss “heavy” particle production processes with the mass ~ \(p\. Hence, our discus¬ 
sion can be also applied to heavy DM production processes. Although in the literature 
it was expected that DM can be produced from the thermal plasma before reheating is 
completed [28, 29], our result implies that they may not be produced in a low reheating 
temperature scenario. Instead, DM may be efficiently produced via scatterings between the 
soft thermalized plasma and hard primaries as studied in Ref. [33], or directly produced 
via inflaton decay. 

In the next section, we briefly explain the finite density effects on boson fields, and de¬ 
rive the relation between the effective temperature which characterizes the finite density 
corrections to the effective potential and the distribution function of radiation. Then we 
investigate the evolution of the thermalization process and calculate the effective temper¬ 
ature of radiation in Sec. 3. In Sec. 4, we apply our results to some important scenarios, 
such as the restoration of PQ symmetry and electroweak symmetry, and the Affleck-Dine 
baryogenesis. Section 5 is devoted to the conclusion. 


1.2 Summary of our results 

Before we explain the detail of our calculations, here we summarize our main results. If 
a system is not in thermal equilibrium, we cannot characterize the system simply by a 
temperature in general. In this paper, we define an effective temperature T* such that it 
describes the finite density effect on the potential of scalar fields because our main purpose 
is to clarify the effects of radiation on scalar fields. Note that these scalar fields of our 

♦2 See also Refs. [31, 32]. 
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interests are different from inflaton, and that inflaton is assumed to reheat the Universe 
via Planck-suppressed decays throughout this paper.*^ See Sec. 2.2 for more details of our 
setup. 

Let us first summarize important equations. If -particles, which are particles inter¬ 
acting with a scalar field (/), are not directly produced from inflaton decay, the thermal 
potential of cf) is roughly given by 






for \g^4) \ < T* 
for \g^(t)\ » r* 


( 1 . 1 ) 


where a^(= g^/A-n) denotes a typical coupling between the scalar field and -fields [See 
Eq. (2.6) and Fig. 1], and a is the relevant coupling constant for thermalization process. 

In Sec. 3, we show that the evolution of the effective temperature for the scalar field 4> 
is dramatically different from that of conventional studies. The evolution of the effective 
temperature T* after inflation is found to be 


- 1/4 
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- ~ < 

m, 


ar 


1/2 


^soft. 


- 1/4 


^4/5p2/5 


for t- • ^ ^ t c 

t-ini ~ ^ ~ ^soft 
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^sort ~ ^ ~ ^max 


for tmax 


( 1 . 2 ) 


where rrij is inflaton mass, a is the relevant coupling constant for thermalization process, 
f j is the decay rate of inflaton normalized by that of dimension-five operator, and t is a 
cosmic time normalized by the inflaton mass. The time scales are given by 


i- ■ = a“^f 

(1.3) 

^soft — ^ ^ 1 ? 

(1.4) 

r — -16/5p-3/5 

^max — ^ ^ I 5 

(1.5) 

f = f-i_E_ 

(1.6) 


where Mpj (^ 2.4 x 10^® GeV) is the reduced Planck scale. We consider the case of < 1 
so that we can neglect non-perturbative effects of inflaton decay. 

To see a rough sketch, we briefly explain the thermalization process after inflation. 
Figs. 7, 8 and 9 maybe helpful to understand the essence of the following discussion. First, 
inflaton decays into hard primaries, whose energy is of order the mass of inflaton m|. At 
the very first stage, a showering of hard primaries produces the soft population. Soon after 
that, the hard primaries then inelastically scatter with each other and release their energy 
into soft particles via in-medium cascading processes. Due to the loss of causality, the 
soft thermal-like distribution with emerges after the time defined by 

If inflaton reheats the Universe via not so small couplings, the preheating and thermal dissipations after 
that may become important. See Refs. [34, 35, 10] for instance. 
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The inelastic scattering rate is suppressed by the LPM effect, which is stronger in lower 
dense environment. Since the density of the hard primaries decreases with time due to 
the expansion of the Universe, the inelastic scattering rate decreases with time. Thus, 
the energy density of soft particles also decreases, so that their effective temperature T* 
decreases with time. Then, at t ~ tgoft, soft particles are completely thermalized. At the 
same time, the number density of the Universe is dominated by the soft particles, so that 
the LPM suppression becomes weaker and weaker. Therefore the effective temperature of 
soft particles increases with time after i ~ Then, at t ~ the primary particles 
can lose their energy completely and are thermalized within the Hubble time scale. After 
that time, the temperature of the Universe can be defined definitely and is determined by 
the energy injected by the inflaton decay.*"^ Since the energy density of inflaton decreases 
with time, the temperature of the Universe decreases, too. Finally, at t ~ inflaton 
completely decays into radiation and reheating is completed. 

There are two things to be noted. 


1. We cannot rely on the conventional estimation of the temperature, T ~ for 

i < tjnax- This is because the energy density is still dominated by the hard primaries 
for i < i^ax’ whose distribution is far from thermal equilibrium. 

2. The maximum temperature of the Universe may not be achieved before ^^d thus 

it may be significantly lower than the standard scenario in which the “instantaneous 
thermalization” is assumed. The effective temperature of the soft sector has two local 
maxima as T^/nii ~ and at the time of t = and t ^^x , respectively. 

In particular, for the case of < a?, the effective temperature reaches its maximal 
value at t ~ t ^ax (see Fig. 11): 

^4/5p2/5^^ for f| < . (1.7) 


If the primary particles was thermalized instantaneously, the temperature of the 
plasma would be calculated from ~ HT^Mpj [28, 29]. (This is shown as 

red dotted lines in Figs. 10 and 11, which clarifies that the thermal effect is overes¬ 
timated in the “instantaneous thermalization” approximation.) The most important 
quantity is the maximal temperature of the Universe because it determines whether 
a symmetry is restored after inflation or not. It is overestimated in the “instantaneous 
thermalization” approximation by the following factor: 


j-Cinst) 
^ max 


I (this work) 
I max 


~ a 


-4/5-p“3/20 



( 1 . 8 ) 


where we assume f| < a^. This implies that the instantaneous thermalization ap¬ 
proximation results in an overestimation by about five order of magnitude for the 
case of a = 1/10, Yj = 10“^^, and Hi = m^, for example. Here note that the decay 
rate of inflaton is bounded from below not to spoil the success of BBN: 




Trh-^pI 
m? 


10 


-27 


V1 MeV J U013 Qevj 


(1.9) 


Still, note that there remains a tail of cascading hard particles as can be seen from Fig. 9. Though this tail 
does not contribute to neither the energy density nor the number density, it can be a source of heavy particle 
production with masses of m » T [33]. 
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Preliminaries 


2.1 Hard primaries 


We focus on the era between the end of inflation and the completion of reheating, dur¬ 
ing which the energy density of the Universe is dominated by that of the inflaton oscilla¬ 
tion. Under the quadratic potential of inflaton, its oscillation amplitude and energy density 
evolve as 


lit) ~ 
Pi ^ 


ait)j 

<ti) \ 

ait)j 


3/2 




( 2 . 1 ) 

( 2 . 2 ) 


respectively. Here, a(t) is the scale factor, is the time at which inflation ends, 1q is the 
initial amplitude of the inflaton oscillation, Tj is inflaton decay rate, and Hj is the Hubble 
parameter at the end of inflation. We define the following dimensionless parameters, f i, 
as: 


K 

i = m,t. 


(2.3) 

(2.4) 


Note that we have i > 1 because the inflaton has to oscillate at least once before it decays. 
We assume that the inflaton decays into radiation perturbatively. This assumption is fulfilled 
at least for < 1. To see this, let us consider the case that the inflaton decays into a fermion 
Xi through the interaction of XlXiXi- This interaction term gives an effective mass to the Xi 
field, so that its frequency is given by = k^+X^I^{t). The non-perturbative decay occurs 
when the adiabatic condition \cb /co^ | ^ 1 is violated. Since i\i^Q < ^ih ^nd k cii m;/2, 
the adiabatic condition can be rewritten as XI q <C rrij. Using the perturbative decay rate 
given by ~ X^rtij, we obtain the condition to the perturbative decay as Tj <C 
Since /g < Mpj, the condition of perturbative decay is always satisfied at least for the case 
of < 1. Otherwise, the inflaton might decay through non-perturbatively [34, 35], which 
is beyond the scope of this paper. 

Since the energy of inflaton decay products (hereafter we call them as hard primaries) 
is of the order of inflaton mass mj, their phase space distribution is given as [36] 

/h(^P) ~ for < p/m^ < 1, (2.5) 

where the factor in the parentheses represents the effect of redshift. It implies smaller 
number density but harder particles than those of thermal distribution for relevant time 
scales in consideration; f^it,mi) ~ <C 1. This is the reason why we call them hard 

primaries. Obviously, number violating processes are essential in thermalization, and in fact 
the hard primaries scatter with each other inelastically and emit soft particles. As we will 
see in the following, these soft particles take crucial roles both in thermalization of hard 
primaries and in dynamics of a scalar field (p which is the main character in this paper. Note 
that the species of soft particles can be different from that of hard primaries. For example, 
even if inflaton decays only into Standard-Model gauge bosons through Planck-suppressed 
operators, soft particles can be quarks and leptons. Thus, almost all of the fields interacting 
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< -► 

Hard Soft 


Figure 1: Schematic figure of our setup: I is the inflaton which creates hard primaries denoted 
by %i, which has typical momentum of mi. Hard primaries emit soft populations of SM particles j 
(including Xi) and also X 4 , (through SM interactions). The field cj) denotes the scalar field that we 
are interested in. It interacts with SM particles via with a typical coupling . 


with a Standard Model particle are affected by the soft sector. In this sense, our calculation 
can be applied to wide range of models. How the soft particles are produced and affect the 
effective potential of (/)-condensate is discussed in the subsequent sections. 


2.2 Setup 

Before going into details, let us clarify our setup. Throughout this paper, we denote the 
scalar field in consideration (which is different from the inflaton) as (p, and fields which 
couple to the scalar field and interact with standard model (SM) particles (e.g., charged 
under SM gauge group) as Xcf,- In addition, inflaton, light fields directly produced via 
inflaton decay, and light fields other than X 4 , are represented by I, Xi, and x, respectively. 
For simplicity, we assume that all the fields except for (p are charged under SM gauge group 
with the predominant coupling being denoted as and that masses of all the particles 
other than X 4 , can be neglected. That is, hard ;fj-particles can emit themselves in general 
and hence x includes Xi, but X(t, is treated separately since it can acquire a large effective 
mass of \g^<p\. Here, symbolically, we represent the coupling between the scalar field <p 
and X 4 , ns g^; concrete examples are*® 

-^Yukawa = g<^4>X^X<PL + h-C., - -^quartic = gl4>^\X4>\^’ <^2.6) 

for instance. Here, x^ii nnd x<p nre fermion and boson, respectively, which are collectively 
denoted as X 4 , hereafter. A schematic figure of our setup is shown in Fig. 1. 

At first, the scalar field condensates almost homogeneously with an expectation value (p. 
Then, when the potential force overcomes the expansion of the Universe, H \/\dV^f(/dcp^\, 
it starts to roll down its effective potential V^ff against the expansion of the Universe. Even¬ 
tually, it relaxes to its “equilibrium value” at that time by dissipating its energy into the 
background plasma. Hence, it is of quite importance to understand the behavior of effec¬ 
tive potential so as to know the cosmological fate of scalar condensates. In particular, the 
effective potential is drastically changed via interactions with a dense and high temperature 
background plasma. 

The field cj) itself can be either charged or singlet under SM gauge group, but note that if (f) is charged 
under SM gauge group, then one should care that SM gauge group is not entirely broken down by the expec¬ 
tation value of (f). Otherwise, the t-channel enhancement of massless gauge boson exchange is suppressed, 
and the following discussion should be altered. 

For a complex scalar field theory, cp is regarded as its radial component. 
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2.3 Thermal potentials 


There are two contributions to the effective potential depending on the expectation value of 
(/)-condensation. The first one comes from the abundant population of -particles, which 
typically has the following form: 


Vm 


sl4>' 




dp“ 


jp 


( 271 )^’ 


(2.7) 

( 2 . 8 ) 


where and (= -^/IpF+ls^^T^) stand for the distribution function and the energy 
of -particles, respectively. For instance, in the case of quartic interaction, g^4>^X^, the 
effective potential for the homogeneous condensate of (p encodes the following term (see 
Fig. 2): 


2 J.2 




p 


) +24(P)] G«(P) gl4>^ ( T [l +24(p)] . (2.9) 


p p 


with = fx^(^p’Py Hsre is the Hadamard/Jordan propagator defined by the 

commutator/anti-commutator: G^^j(x,y) = Tr / [j^(x), ;f_j^(y)] V In the first similarity. 


we have assumed the Kadanoff-Baym ansatz between the Hadamard and Jordan propaga¬ 
tor, and then assumed that the spectrum is dominated by particle-like excitations in the 
second similarity. One can see that the finite density correction is reproduced aside from 
the quadratic divergence that should be canceled by the counter term. A similar computa¬ 
tion yields the finite density correction from the Yukawa interaction, g(j,4> X 4 >rX 4 >L’ which is 
essentially the same as Eq. (2.7). The expression ofEq. (2.7) coincides with the well-known 
thermal mass, if is given by the thermal distribution and if the effective mass of x,p is 
smaller than the temperature, ~ |g^(/)| <C T. Motivated by this observation, we define 
the effective temperature for -particles as follows: 



r fx^p^ 



( 2 . 10 ) 


Although the distribution might be far from thermal equilibrium, we can estimate the 
effective potential for the scalar field (p, which is caused by abundant -particles. Imitat¬ 
ing the expression at thermal equilibrium, we write down this effective mass contribution 
as 


Vm 





( 2 . 11 ) 


where = gj/47i. 

When the field value of (/)-condensate is larger than the typical production scale of 
-particles, they are less likely to be produced from other light particles. Hence, the 
contribution given in Eq. (2.7) maybe suppressed. Even in this case, the effective potential 
for the 0-condensate receives finite density corrections from other light particles. This 
second contribution can be understood from the fact: the runnings of coupling constants are 
affected by the expectation value of the scalar field. Suppose that (p has a large expectation 
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Figure 2: Diagram responsible for the finite density correction to the mass term of cj) in the case of 
quartic interaction. 


value. Since Xcj, has a large effective mass of\g^(p \ and is less likely to be produced, one can 
integrate out -fields. As a result, the running coupling constant of gauge group under 
which X 4 ) is charged encodes the (/)-dependence as a threshold correction of -fields, 
or in other words, the (/)-condensate interacts with other light particles via the operator 
~ Here we have assumed that cf) is singlet under this gauge group. 

This implies the following contribution to the effective potential for the (/)-condensate (see 
Fig. 3):*® 





g\4>)T:. 


( 2 . 12 ) 


Here we have defied the effective temperature F* for light particles x, which does not 
always coincides with that for -particles since they have masses of |g^0|. Extracting 
0-dependent part, one may obtain 







(2.13) 


with a = g^/A-n. This is nothing but the so-called thermal-log potential for the thermal 
distribution of -particles [27]. 

Thus, in order to discuss the finite density corrections to the effective potential for (j), 
we have to know distribution functions for light -particles as well as that for -particles 
with the mass oi\g^(j)\. In the following, we discuss how x and Xcj) are produced from hard 
primaries during the course of thermalization and calculate their effective temperatures T* 

Before closing this section, we would like to comment on basic assumptions behind 
our following discussions. First of all, the effective potential listed above is computed per- 
turbatively which implicitly assumes that the occupation number /(kmax) [that dominates 
the integral /^//k ~ /(kmax)k^ax] should be smaller than a~^, since otherwise the per¬ 
turbative expansion breaks down. Below, we will see that this condition is automatically 

Other interactions like Yukawa may also encode the threshold correction of . 

Here we have assumed that is also charged under the dominant gauge group of coupling g. If this is 
not the case, a factor should be multiplied, which is defined as the ratio of the coupling under which is 
charged to the dominant coupling g. We restrict ourselves to order of magnitude estimation and do not care 
about this factor in the following discussion to avoid complications. 


9 





Figure 3: Diagram responsible for the logarithmic correction to the effective potential of cj). 


satisfied. In addition, although we only focus on the effective potential to discuss the sym¬ 
metry restoration, it is not trivial that whether or not the effective potential is maintained 
after the onset of oscillation and also the scalar condensate soon dissipates its energy into 
background plasma. In the following, we simply assume that cf) can dissipate its energy 
into background plasma soon after the onset of its oscillation, since its relaxation process 
can be strongly model dependent [37]. 


3 Thermalization History 

In this section, we investigate the thermalization process of hard primaries and calculate 
the effective temperatures for the soft sector T* and T* . Basic parameters which are re¬ 
quired to study thermalization after inflation are the mass of inflaton m^, the decay rate of 
inflaton and the dominant couplings of decay produces a; e.g., for the SM plasma, the 
predominant interaction is the strong coupling, a = a^. This is because the first two param¬ 
eters nil, f / determine the typical distribution of decay products, and the last parameter a 
is essential to estimate their interaction time scale. We are also interested in the production 
of -particles and hence the coupling between Xtf, ^nd the scalar condensate (j), which is 
given by g^, is also an important parameter. Thus, we have four parameters essentially to 
capture this system. 

We assume that the inflaton decays into particles Xi which are charged under the abelian 
and/or non-abelian gauge groups (SM gauge group). We mainly consider the era after 
i = tjnj, where is defined below [see Eq. (3.22)]. First, the hard primaries emit soft 
particles but the soft sector is not thermalized soon. As shown in Sec. 3.2, the soft sector is 
completely thermalized after the time of t = [Eq. (3.43)] . After that, around the time of 
i = [Eq- (3.57)], the hard sector as well as the soft sector are completely thermalized 
and the effective temperature reaches a (local) maximum. Then the reheating is completed 
at t = [Eq. (3.63)]. We investigate thermalization processes during these eras in turn, 
and then we briefly summarize the result of this section in Sec. 3.5. Our procedure is based 
on the one used in Ref. [38, 39] (see also Refs. [30, 33]). 

Figs. 7, 8 and 9 summarize the essence of this section. Hurry readers can skip all the 
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details and move to these figures. 


Kinetic equations 

Here we briefiy summarize minimal knowledge required for our estimation. Eventually, we 
will see that the collinear splitting process in medium plays crucial roles in the following 
discussion. We also summarize basic formulae of this process in SU(N) with Np-Ravor 
fermions in Sec. A for the sake of completeness. 

Our method is based on kinetic equation (Boltzmann-like equation), which can be de¬ 
rived from first principles, i.e., from Schwinger Dyson equations in the closed time path 
formalism (see for instance [40, 41], reviews [42, 43] and references therein). It was 
shown that particle-like excitations of mj <C p with being the screening mass obey the 
kinetic equations if these quasi-particles fulfill the following conditions: (i) Typical “size” 
of quasi-particles are smaller than the mean free path, (ii) Typical duration time of interac¬ 
tions are shorter than the mean free time, (hi) Typical value of distribution function should 
be smaller than 1/a. The condition (ii) is sometimes violated for collinear emissions in 
the medium, and in this case appropriate resummations are required, which is so-called 
the Landau-Pomeranchuk-Migdal (LPM) effect [44, 45, 46, 47, 48, 49, 50]. After the re¬ 
summations, we eventually obtain the following kinetic equations for quasi-particles in the 
plasma [51, 38, 39]: 


d^—Hk 


dk 




(3.1) 


where /, is the distribution function of quasi particles with • being species and is the 
collision term. As we will explain later, at least for the perturbative Planck-suppressed decay 
of infiaton, these conditions are safely satisfied*^ except for the condition (ii). Hence, we 
can rely on the kinetic equations given in [51, 38, 39] . In our case, we are interested in the 
distributions of x and x<i) to obtain information of effective temperatures. 

The collision term for -particles encodes two-to-two scatterings responsible for diffu¬ 
sions and effective “one-to-two” (inelastic) scatterings due to splittings: 

^[/J=>(^2to2 + ^split)[/J, (3.2) 


where 


2 to 2 






b,c,d -J 


p,k',p‘ 


, 2E^.2Ep2E^.i2E„i 


(271)^5^^^ {K + P-K' -P') 


X [/a(fc)/b(p) [1 ±/c(fcO] [1 ±/d(pO] - (inverse process)] , 

(3.3) 


It is contrary to the preheating case where the fluctuations can be strongly correlated. See Refs. [52, 53] 
for instance. 
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^split [/a] 


(2nf 


^ dp'dk'5 {k-p' - k') yl, {k; p', k') 


b,c -J 


X [faW [1 ±/b(pO] [1 ±/c(fcO] - (inverse process)] 
3 r 


+ <ip'<ip5ik + p-p')rl^{p';k,p) 

^ b,c j 

X [faWfbip) [l±/c(pO] - (inverse process)] 


(3.4) 


dlnp'r,piit(pO [4(fc) [1 ±/.(fc-pO] [1 ±/.(p 0] - (inverse process)] 


+ 


dlnp' 1^—J r,piit(fc) [/;^(fc)/.(p^ - k) [1 ±/.(p 0] - (inverse process)] , 

(3.5) 


with Vq being the number of degree of freedom for a species a (normalized by one real 
d.o.f.). Here, the signs are taken as + (—) for boson fields (fermion fields). FspUf stands for 
the splitting rate. One can show that Eq. (3.5) can be derived by using the definition of the 
splitting function, for a be, given in App. A. Instead, below, we give an intuitive 
derivation of the splitting rate, figpiit, in order to understand its physics. See App. A for more 
rigorous explanation. 

The first term of the right-hand side of Eq. (3.2) imprints elastic scatterings that diffuse 
distribution functions in momentum space, where the process is dominated by small mo¬ 
mentum exchange due to the strong t-channel enhancement of gauge bosons as shown in 
Eig. 4. The momentum transfer squared obeys the following diffusion equation: 

(Ap)^~qeit, (3.6) 


with 


Pel 


, 


dy ci V • 


/.(pO[i±/.(pO]. 


(3.7) 


Here we have used the rate of elastic scatterings 


dr 


el 


a 


dqi qiiqi + m^) 


/.(pO [i±/.(pO]. 


(3.8) 


Note that the screening mass for -particles is given by ~ clT^. 

Then, let us move to the second term, ‘T^gput. As notified in the beginning, in estimating 
the inelastic scattering rate of the hard primaries, it is necessary to take into account an 
interference effect between a hard primary and its daughter particle, known as the LPM 
effect. A Eeynman diagram describing the interference effect is shown in Pig. 5. The inter¬ 
ference effect forbids subsequent scattering processes during the interval while their phase 
spaces overlap with each other. Taking this effect into account, we obtain the following rate 
of inelastic scatterings: 


rspiit(fc) ~ « Min 



1 ■ 

^form _ 


(3.9) 
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Figure 4: Diagram for elastic 2 to 2 scattering process. 



Figure 5: Diagram representing an interference effect, known as the LPM effect. 


with 


tform='l/—• (3-10) 

V *3el 

Here we omit a factor that comes from model dependence and logarithmic enhancement. 
See App. A. The time scale of denotes the formation time which it takes to resolve 
the overlap between the daughter and parent particles, k is the momentum for a daughter 
particle. One can see that the rate is LPM suppressed above the threshold momentum /clpm, 
which is given by 


^LPM — p2 • (3.11) 

Note that Eq. (3.10) comes from the condition so that the overlap between the parent and 
daughter should be resolved; t > k/k^ with k^ ~ q^it. Therefore, for a given time t, there 
is an upper bound on the daughter momentum: 

fcform = (3-12) 

Above fcformj even the diffusion cannot resolve the destructive interference. In this case, the 
vacuum cascades may become important, though the emitting angle between the parent 
and the daughter is bounded below 6 > {kt)~^^^ for a given k and t. See also Eq. (3.28) 
and discussion about it. 

So as to obtain the effective temperature of X 4 ,, we also have to know how -particles 
are produced. The collision term for -particles responsible for -production maybe given 

by 


^ (“^210 2 + “^split + ^dec) (3.13) 

where the last term with the subscript “dec” denotes the decay/inverse decay processes. 
They are expected in general since the 0-condensation can give a sizable mass to x<i,- 













X(l> 

-► 

Figure 6: Diagram for the decay of /^-particles. 



particles; that is, larger than the screening mass, \g^4>\ > m^, and are calculated by 


“^dec \-fx^ ] 


decv 


k',p' 2£fc2£(./2£p/ 


(271^5^4^ {K-K' -P') 


[1 ±/.(£')] [1 ±/.(p0] - (inverse process)] . (3.14) 




These collision terms depend on models how -particles couple with other light particles; 
in particular, we have to specify interaction terms which can induce decay of Xcj, into light 
particles for a sizable mass of Xcf), \s<i>4’\ > nr^. For simplicity, we assume that the typi¬ 
cal magnitude of this term is proportional to e^a with e being a small parameter and a 
being the fine structure constant of the gauge group which dominates the thermalization. 
We assume < a to avoid unnecessary complications. See also discussion given below 
Eq. (3.34). For instance, if -particles are charged under the gauge group and if they mix 
with light particles, then their decay rate may be proportional to e^a\g^ (p\, where e is iden¬ 
tified as the mixing angle (see Fig. 6). Under this assumption, the dominant contribution 
to first two terms is essentially the same as that for -particles, though we should note that 
when a daughter particle has a mass larger than the screening mass m^, we have to replace 
rUj in the scattering rate of Eq. (3.8) to the mass of the daughter particle. And also note 
that these terms can be sources of -particle production. 


^ini ^soft 

Apparently, two to three processes with the momentum exchange of the order of rtij seem to 
be efficient in order to decrease/increase energy/number of hard primaries. However, this 
is not true since the reaction energy/occupation number is too large/small and the rate is 
found to be much smaller than the Hubble parameter. This is a typical situation of reheating 
via Planck suppressed decay, < 1. In this case, the hard primaries lose their energy 
mainly through inelastic collinear scatterings with the rate of Eq. (3.9), which are enhanced 
by the t-channel contribution [31]. Below, we will see how -particles are produced via 
inelastic scatterings, and form the soft population. Then, we discuss production of X(i>- 
particles, which is essentially different from that of -particles due to the mass of Ig^^l. 


Production of -particles 

At the first stage of thermalization process, the hard primaries scatter among themselves 
and emit soft particles through the inelastic scatterings imprinted in ‘^sput- The phase space 
distribution of these soft particles might be estimated as 

/,(t,k) ~ r,piit(k)n;,k“^t, (3.15) 
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where ~ J is the number density of hard primaries [see Eq. (2.5)]. However, low- 
momentum modes are so over-occupied that inverse processes and subsequent elastic scat¬ 
tering processes have to be taken into account. In fact, soft particles feel subsequent scatter¬ 
ings with hard primaries*^® and low-momentum modes fall into a thermal-like distribution 
such as*^^ 

~ ^ for k < k^iax. (3-16) 

k 

where is determined by the diffusion constant of Eq. (3.7) with f = fh and is written 
as 

fcmax~ (3.17) 

One can show that it is the same as the threshold momentum of LPM suppression, i.e., 
fcmax ~ ^LPM- T^e effective temperature for bosonic soft modes can be estimated from the 
energy conservation. The energy of soft modes below the scale k^^^ is given as*^^ 

fcmaxrsplit(fcmax)tflh ~ Ps ~ ~ (3.18) 

where E^piij is given by Eq. (3.9). Thus, the effective temperature for the soft modes is 
given as 

T, ~ (3.19) 


from Eq. (3.18) together with Eqs. (3.9) and (3.17). Here, we have used TspUt ~ oi\f^^Jk 
because the relevant energy scale k^^^ is so large that such a soft mode are produced by 
LPM suppressed interactions. 

We illustrate the distribution of -particles in Pig. 7, where we include the hard sector 
as well as the soft one. The distribution of smaller wavenumber modes than k^^ is a 
thermal-like one given by Eq. (3.16), while that of larger wavenumber modes in the soft 
sector is determined by the LPM effect as Eq. (3.15). The distribution of hard sector is 
given by Eq. (2.5). The distributions of soft and hard sectors coincide with each other at a 
certain wavenumber, shown in Pig. 7. Por clarity, we call the distribution of 

-particles above/below this momentum as the hard/soft sector. 

At this regime, the effective temperature T* is dominated by the contribution from the 
soft sector 


t -'2 
i , ~ 


> 


fs 


r^Tk ~ 

^s^max ^ ^ I ^ 




(3.20) 

(3.21) 


Although there are many processes for the thermalization of soft modes (see Ref. [38, 39]), we explain 
one of them as an illustration. 

In the case that soft particles are fermions, the distribution A cannot exceed unity due to the Pauli 
blocking effect. In this case, we should truncate when it reaches unity. However, we expect that soft 
particles contain bosons that dominate the thermal effect on c^, so that we focus on the distribution of bosonic 
particles. Its extension to fermionic case is straightforward. 

One might wonder if the energy conservation within the soft modes [that is, Eq. (3.18)] is violated due to 
the scatterings between hard primaries and soft modes. However, this does not change our order estimations. 
The energy for soft modes is dominated by that for the modes around fcniaxj ^^d the scatterings between hard 
primaries and such soft modes are inefficient to change the energy-conservation relation of Eq. (3.18). 
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• ~ ~ ~ ~ 1 /2 
Figure 7: Distribution of j-particles for tj^j < t < Recall that fcj^ax ttt^, fcform ~ 

fcmax(t/tini) tjni = a“ ' . Soft (pink): The thermal-like distribution dominates below 
the LPM-suppressed spectrum can be seen between kj^ax and fcform> and the vacuum cascade may be 
relevant above kform- Hard (blue): The hard spectrum is sourced by the direct decay of inflaton and 
by its red-shifted spectrum of previously produced one. See Eq. (2.5). 


where the inequality holds for i > (see below). This implies that the 

screening mass (= aT^) is already dominated by the soft sector. Note that ^ T* 
because the soft sector is not completely thermalized and and T* are nothing but effective 
temperatures with different definitions. 

Here, we should ensure that the relevant wavelength scale (e.g., is smaller than 

the Hubble horizon ~ t for a consistent treatment because longer wavelength modes 
than the Hubble horizon cannot be emerged due to the loss of causality. This implies that we 
can determine the time scale when the softest mode emerges by comparing the maximum 
momentum of the soft sector with the Hubble parameter. The inequality > H indicates 

i > (3.22) 


This inequality immediately means that the occupation number of soft sector at k^^^, which 
dominates the energy and number density of the soft sector, is smaller than a~^: 


f(tk ) ~ 

''‘max/ 


a 



(3.23) 


Hence, we may treat the soft sector perturbatively. Then, let us compare the screening mass 
with the maximum momentum kmax- Up to here, our arguments rely on the effective kinetic 
theory of non-Abelian plasma, and hence produced modes should not be too soft, k > m^, 
for a consistent treatment. In fact, Eq. (3.22) also ensures 


We should also note that 



(3.24) 


(3.25) 
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which means that the screening length is smaller than the horizon length. Finally, we check 
that a typical time scale of interactions (e.g., r“j^) is smaller than the cosmological time 
scale 






(3.26) 


where we use Eq. (3.8). Therefore, the condition t > ensures that > H, 

/j < a~^, and Tgi > H. These inequalities justify our calculations performed in this paper. 
Otherwise we may have to take into account the effect of causality, the finite cosmological 
time scale, and non-perturbative effects on production processes. 

Finally, we comment on the era before tj^i. Since elastic scatterings among decay prod¬ 
ucts does not take place within the Hubble time scale, medium induced cascades can be 
neglected. In addition, the formation momentum, kformj is smaller than the Hubble param¬ 
eter before tin;: 


^form 



(3.27) 


Still, the vacuum cascades might broaden the spectrum towards the infrared, and hence 
we roughly estimate their effect on the effective temperature. The quantum formation time 
scale in vacuum implies that one can find a quanta of k with a probability of &{a) per 
lnkln0 for 6 > {kt)~^^^. Omitting the log factor, one may obtain the distribution of the 
soft sector and the effective temperature as 


/,(t,k)~af^t ^ 



for k> H 





(3.28) 


which coincides with Eq. (3.20) for i —>• Hence, we expect that the effective tempera¬ 
ture is at most that of After these particles with k~^ are soon red-shifted away and 
LPM suppressed spectrum starts to dominate the soft sector. Thus, we will not consider the 
case of t < ti„j further in the following. 


Production of -particles 

In contrast to -particles, -particles can have sizable masses of due to the large 

expectation value of the 0-condensate. The distribution in momentum space can be dif¬ 
ferent from that of -particles and hence separate discussions are required. Although the 
calculations are complicated, we find that the effect on the effective potential from x 4 ,- 
particles is always subdominant for the case oi \g^(j)\ > a~^k^^ and f| < a®. Even in 
the case of f ^ > a^, its effect is roughly the same order with that of thermal log potential 
from -particles, so that we can neglect it for rough estimations. We briefly explain X 4 - 
particle production processes in this subsection, while detailed calculations are shown in 
Appendix.*^^ 

♦13 yjq conservatively omit the decay of massive -particles to show that the effect on the effective potential 
from X 4 -particles is usually subdominant. In other words, we show that the effect from X 4 -particles can be 
neglected even if they are stable and maximally abundant. Note that the calculation can be applied to the 
production of massive dark matter, once we forget the field 0 and identify X 4 as just a dark matter with mass 
of |g^0|. Its abundance can be calculated from the resulting distribution function of -particles given in 
Appendix. 


17 





If the effective mass, \g^(p\, is smaller than the screening mass, m^, then the results are 
exactly the same as those for -particles. Hence, we concentrate on |g^0| > in the 
following. First, let us consider contributions from splittings which are imprinted in 
Note that since we consider the case of |g^0| > m^, we should replace in Eq. (3.8) 
to \g(j,<P\- The phase space distribution and the resultant effective temperature may be 
estimated as 


where 




hard 


rspiit(fc)fthfc 


-3. 


'■,X4, 


hard 




(3.29) 



(3.30) 


The “hard” subscript means that splittings occur dominantly through interactions among 
hard primaries, since the number density is still dominated by the hard sector at this stage, 
i < tjoff. Here Tsp^j is given by Eq. (3.9) but is replaced by |g^(/>|. In this case, the 
threshold momentum is given by 


^LPM,x^ “ ^LPM 



(3.31) 


If the effective mass for X 4 , is not so large, that is, (g^(/>)^ < then the distribution 

function for X 4 > coincides with that for x below . 

In addition, -particles are produced via ‘ 1^2102 with interactions among the 

soft sector and also between the soft and hard sector. Eirst, let us concentrate on "^2 to 25 
since it gives model independent contribution.*^"^ Although the number density is still 
dominated by the hard sector, the soft sector may produce “soft” -particles with a larger 
rate (e.g., via s-channel scattering processes). The effective temperature may be given by 



soft 



dlogic' nXk')n, 



(3.32) 


for \g^4)\ > where n, ~ j^f,. Note that, for \g^(f)\ < k^^^, it is given by 
k^ . They are also produced from interactions between the hard and soft sectors: 


'X 4 , 


soft 


s/0 


a 


\g6(t>\' 


dlogic' nXk')n^ 


\gs4>\^ 


(3.33) 


where the subscript “int” indicates contributions from the interaction between the hard 
and soft sectors. Note that both and can be calculated from the distribution function 
illustrated in Eig. 7. Then, we give comments on contributions from which strongly 
depends on details of a model. As an illustration, we have assumed that the magnitude of 
this process is proportional to e^a. The inverse decay can produce ;f^-particles and it may 
yield 



soft/int 





dlogk\(k'X/^ 



♦14 (jj-Qp terms proportional to at least for < a. 


(3.34) 
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Hence, if the typical mixing parameter e is not so large, say < a, then this contribution 
does not exceed those from Eqs. (3.32) and (3.33). Moreover, this collision term, 
yields decay of -particles: 



~ e^a\g^4)\t < 




(3.35) 


where we use defined in the next subsection. One can see that the decay 

is insignificant for \g^<p\ < /Cj^ax at this stage, i < tso^. As explained in the footnote 4^13, 
we can omit the decay term in order to show that the effect on the effective potential from 
abundant -particles is subdominant for |g^0| > 

For later convenience, we denote all these contributions where -particles are pro¬ 
duced indirectly (not directly produced by inflaton decay) as 


= max 


, T ^ 



indir 

*,X4, 

hard 

soft 

int. 


The explicit forms of the RHS are shown in Appendix. 


(3.36) 


In contrast to the above discussion, if -particles are produced directly from inflaton 
decay, then the situation turns out to be slightly different. This is because the hard sector 
still has a sizable contribution to the effective temperature, which tend to dominate for a 
large field value. In this case, its distribution function has a contribution of Eq. (2.5). As 
shown below, the effect on the effective potential from -particles can be efficient com¬ 
pared with the thermal log potential. Therefore, in this case, we should take into account 
the decay of -particles to obtain more realistic predictions because -particles obtain 
effective masses of \g^(j)\ and can decay into light particles in general (see footnote 4^13). 
Its decay effect can be taken into account by multiplying Min[l, (T^jg^-ayt)"^] to its distribu¬ 
tion of Eq. (2.5), where T^gj-gy is the decay rate of -particles. The effective temperature 
may be estimated as 


where E^ecay ~ 



dir 


= f^t ^m^Min 



(3.37) 


Thermal potential 

Here we summarize the thermal potentials given in Eqs. (2.11) and (2.13), and discuss 
which contribution dominates for each case. See also Appendix B. Note that for an interval 
between < a^^kj^ax, the effective temperature may depend on e^. We do not 

care about this small interval for our rough estimation. 

• Indirect -production (produced from hard primaries/soft daughters): 

- The simplest case is a large field value regime, > a~^k max - In the case 

of fj < a^, the effective mass from abundant -particles always gives sub¬ 
dominant corrections to the effective potential, meanwhile the log contribution 
dominates: 
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Even in the case of f ^ > a^, the effect of -particles is very limited, so that we 
can neglect its contribution for a rough estimation. 

- For a small field value regime, < k^ax, the effective mass from abundant 

-particles dominates the effective potential: 


■,X4, 


indir 


5a4,T^; 


5 e [a, 1]. 


(3.39) 


Here we do not write down complicated parameter dependence of T^ ^^^lindir 
brevity. See Appendix B for details. 


• Direct -production (produced from inflaton decay): 

- For the large field value regime, \g^4> \ > the effective potential is given 

by the competition of two contributions; the log dependence of the running 
coupling constant and hard -primaries from the direct decay of inflaton. The 
contribution from the direct decay of inflaton, which is given by 



dir 






1 


r t 

^ decay ^ 


dominates the effective temperature for 


ant;Max 






(3.40) 


(3.41) 


- For the small field value regime, \g^(p\ < k^^, the contribution from the abun¬ 
dant -particles dominates over that from the running coupling constant, as 
explained in the previous case. Again, we do not write down the explicit pa¬ 
rameter dependence of T*_;^^lindir since it is complicated. See Appendix B for 
details. 

The trivial case is for there the contribution from abundant X(i>- 

particles always dominate over that from hard X 4 , '■ 


a. r 




indir 




X4, 


dir 


(3.42) 


This is simply because the effective temperature for massless modes is domi¬ 
nated by the soft sector as shown in Eq. (3.21). 


^soft ^max 

After the time when the effective temperature T* becomes comparable to k^^^, the soft 
sector is completely thermalized by their own interactions and can be described by a single 
parameter Tj ~ T* ~ ^max- This occurs at the time of 

= (3.43) 

At the same time, the number density of soft sector is as large as that of hard primaries, 
which means that we should substitute / in Eq. (3.7) with f^. 
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Production of -particles 

The energy injection to the soft sector is given by 

Ps ~ Kpik^h- (3-44) 

The momentum scale is given by a criteria r(fcspii() ~ H(t), which means that modes 
below fcgpiit are soon thermalized and fall into soft modes. Equating with T^, we obtain 

T.r-. (3.45) 

fcspiit ~ (3.46) 

where we use ~ The thermal potential is now determined by the temperature of 

the soft sector: 

^max ~ T* ~ T, ~ (3.47) 


Note that the distribution function for soft -particles after i > is given by 


f 


k 


for k < T„ 


fsit,k) ~ 



for T,<k< fcform, 


af|t ^ 

I 



for 


(3.48) 


Also note that the “hard” interactions among soft particles with the momentum exchange 
of is faster than the cosmic expansion: 


H 



> 1 . 


(3.49) 


Hence, the soft sector is thermalized separately. 

We illustrate the distribution of -particles in Fig. 8, where we include the hard sector 
as well as the soft one. The distribution of the soft sector is given by Eq. (3.48), while that 
of the hard sector is given by Eq. (2.5). Eventually, the distributions of soft and hard sectors 
coincide with each other at the wavenumber of fc ~ where 

is defined in the next subsection [See Eq. (3.57)]. 


Production of -particles 


If the effective mass term, \g^(p\, is smaller than the temperature of the soft sector, \g^(p\ < 
Tj, then the dominant contribution to the distribution function of -particles is given by 

fx^it,k) ~ for fc < T„ (3.50) 

which yields 



indir 




(3.51) 
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Figure 8: Distribution of -particles for tjoft < t < Recall that T* rs^ ar/ m^(t/tsoft), 

fcform ~ «“^7’*(t/tsoft)^, tsoft = and t^^ax = [See Eq. (3.57)]. Soft (pink): 

The thermal-like distribution dominates below r*(= Tj), the LPM-suppressed spectrum can be seen 
between ^nd fcform> ^nd the vacuum cascade may be relevant above fcform- Hard (blue): The 
hard spectrum is sourced by the direct decay of inflaton and by its red-shifted spectrum of previously 
produced one. Eventually, the distributions of soft and hard sectors become to coincide with each 
other at the wavenumber of 


for Ig^^l < Tj. Basically, this is because “hard” interactions among the soft sector with the 
momentum exchange of is much faster than the cosmic expansion at this stage, as shown 
in Eq. (3.49). Thus, -particles with a mass lighter than T. are efficiently produced and 
participate in the thermal plasma. 

The case of |g^(/)| > T^ is calculated in Appendix. The -particle production processes 
are the same with the ones in the previous subsection, though the number densities of 
and are different. 

Thermal potential 

Here we summarize the thermal potentials and discuss which contribution dominates for 
each case. See also Appendix B. Again, note that for an interval between T* < \g^4>\ < 
a~^T^, the effective temperature may depend on e^. We do not care about this small interval 
for our rough estimation. 

• Indirect -production (produced from hard primaries/soft daughters): 

- First, we consider a large field value regime, \g^4>\ > In the case of < 

a^, the effective mass from abundant -particles always gives subdominant 
corrections to the effective potential, meanwhile the log contribution dominates: 

Even in the case of f ^ > a^, the effect of -particles is very limited, so that we 
can neglect its contribution for a rough estimation. For i > Min 
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that is, soon after the soft sector is thermalized the running coupling constant 
contribution dominates for |g^0| > a~^T^ independently of fj. This is because 
the number density of soft sector increases with time, while the number density 
of hard primaries decreases due to the red-shift of inflaton number density. 

- For a small field value regime, < T*, the effective mass from abundant 

-particles dominates the effective potential: 






indir 


(3.53) 


• Direct -production (produced from inflaton decay): 

- For the large field value regime, the effective potential is given 

by the competition of two contributions; the log dependence of the running 
coupling constant and hard -primaries from the direct decay of inflaton. The 
contribution from the direct decay of inflaton, which is given by 


'■,X4, 


Min 




'-soft 


s<j> 4^ 


(3.54) 


dominates the effective temperature for 

/ r \ 5/2 


rrij > 


> am 


I 


Min 


^soft 


a-f-a - 


^soft J 


7/2' 


(3.55) 


- For the small field value regime, \g^(p \ < T*, the contribution from abundant x,p- 
particles is always larger than other ones; from the running coupling constant 
and also form the direct decay of inflaton: 


V *,X4, v *,X4, 




indir 


*,X 


(3.56) 


This is simply because the effective temperature for massless modes is domi¬ 
nated by the soft sector as shown in Eq. (3.21). 


3*3 tmax < t < tRH 

When r 3 piit(fc ~ m^) H is satisfied, that is, when becomes comparable to m^, primary 
particles lose their energy completely and are thermalized soon. This occurs at the time of 

= (3.57) 

This is the actual thermalization time scale of hard particles. 

Production of -particles 

The energy conservation implies that 

p,~p^r^t~H(t)r^M2j, (3.58) 
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Figure 9: Distribution of -particles for t ^av < t < tR^. Recall that T* ~ 

tmax == , and tj^ = ^9- (3.63)]. Note that there remains the tail 

of LPM-suppressed spectrum oc because the inflaton continuously produces primaries with the 

momentum of p ~ m^. 

where (~ is the energy density of radiation. This gives us the time evolution of the 
temperature such as 

r* ~ (3.59) 

Since the hard sector as well as the soft sector are thermalized, they can be described by a 
single parameter T*. 

We illustrate the distribution of -particles in Fig. 9. The distribution of smaller wavenum¬ 
ber modes than T* is a thermal one, while that of larger wavenumber modes is determined 
by the LPM effect as Eq. (3.15). Note that the hard sector is completely thermalized, so 
that it is absent in Fig. 9. 


Production of -particles and thermal potential 


The production process of -particles is the same with the one considered in the previous 
subsection. However, since the number density of hard primaries decreases efficiently by 
splittings for t > the contribution from the direct decay of inflaton is now given by 




~ frt ^ Min 


dir 




(3.60) 


Next, we discuss which contribution dominates the thermal potential. See also Ap¬ 
pendix B. Again, note that for an interval between the effective 

temperature may depend on e^. We do not care about this small interval for our rough 
estimation. 


• Indirect -production (produced from hard primaries/soft daughters): 

- For the large field value regime, \g^4>\ > we find that the effect on the 

effective potential from -particles is always subdominant compared with the 
thermal log potential of Eq. (2.13). 
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- For the small field value regime, < T^, the distribution of -particles and 

its effective temperature are given by Eqs. (3.50) and (3.51). In this case, the 
effective potential of <p is given by the usual thermal mass term of Eq. (2.11). 


• Direct -production (produced from infiaton decay): 

- Eor the large field value regime, the effective potential is deter¬ 

mined by the competition between the contribution from the running coupling 
constant and the direct decay of infiaton. The latter yields 




dir 


(afj'q 


8/5 


9/8 


r.^Min 




1 / 8 - 


This contribution dominates the effective potential for 


5/16 


m, > 


> am 


I 


Max 


1. 


1/2"\-6/5 




11 / 16 ' 


t J 

(3.61) 


(3.62) 


- Eor the small field value regime, < T^, the effective potential of cj) is 

governed by the usual thermal mass term of Eq. (2.11). 


3-4 Irh < t 


When Tj H{t) is satisfied, the energy density of radiation becomes to dominate that of 
the Universe. This is when the reheating is completed: 

t2 

(3.63) 


r r _ f5-l pi 
f ~ ^RH — 2 ■ 


The reheating temperature is given by 

T’rh ~ (3.64) 

After the reheating is completed, the temperature of the Universe decreases with time as 


1/2 


T* ~ ^’rh 


'-RH 


(3.65) 


because in the radiation dominated era. 


3.5 Summary 


To sum up, we obtain the evolution of the effective temperature of -particles as follows: 
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(3.66) 
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Figure 10: Evolution plots for the effective temperature of the soft sector T*. We assume a = 0.1, 
Hj = rrij, and rrii = 10^^ GeV. We take fj = 1, which corresponds to the case of fj > a^. The blue 
lines are results derived in this paper, while the red dotted lines describe the temperature derived in 
the literature by assuming the “instantaneous thermalization”. We also show the effective tempera¬ 
ture of the hard sector as green dashed lines, for the case that they are produced from inflaton decay, 
where we neglect the decay of -particles. The blue dot-dashed line is the effective temperature 
for i < tjjii, which is estimated from the discussion of the vacuum cascades [see Eq. (3.28)]. 

where the time scales are given by 

(3.67) 

(3.68) 

(3.69) 

(3.70) 

Here note that we have implicitly assumed tj^ > ^max so far, which implies the thermal- 
ization takes place faster than the complete decay of inflaton. As shown in Ref. [30], this 
inequality is satisfied in most cases: a > 4 x 10“^(mj/10^^ GeV)^^®fy®. These results are 
illustrated in Figs. 10 and 11 as blue lines. The effective temperature of the soft sector 
has two local maxima as T^/rUj ~ and at the time of t = tj^j and 

respectively. The global maximum of the effective temperature is given by 


ai/^fy^mj for 1 > fj 
^4/5p2/5^^ for fj < 


(3.71) 


There is a local minimum given as T* r\^ afat i ~ tgoft. Note that the effective tempera¬ 
ture for i < tjni is at most that for i ~ tj^i from the discussion of the vacuum cascades [see 
Eq. (3.28)], so that the maximal temperature is in fact given by Eq. (3.71). The thermal- 
ization process described in this section is summarized in the second paragraph in Sec. i.2. 
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Figure 11: Same as Fig. 10, but we take f/ = 10 which corresponds to the case of f/ < a^. 


Now we obtain the evolution of the thermal effects on the potential of (/). If -particles 
are not directly produced from inflaton decay, the thermal potential of (p is roughly given 

by 






for \g^(p\< r* 
for » r*. 


(3.72) 


for the case of < a®. Here, we neglect the complicated parameter dependence of the 
effective temperature for aT* < \g^4>\ and roughly evaluate it as ~ T* 

for \g(j)4>\ <C T*. Even in the case of f| > a®, the deviation of the above formula is very 
limited, so that one can still use it for a rough estimation. On the other hand, if X(i>- 
particles are directly produced from inflaton decay, it may affect the thermal potential of 
(p. In particular, for the large field value regime, \g^<p\ ^ the effective potential is 

given by the competition of two contributions; the thermal log potential from -particles 
and thermal mass from -particles. The latter contribution is given in Eqs. (3.40), (3.54) 
and (3.61). If we can neglect the decay of -particles owing to the smallness of e, 
evolves as green dashed lines in Eigs. 10 and 11. 

Einally we comment on the uncertainty of our results. We omit the log-enhancement 
factor of splitting rate in Eq. (3.9) for simplicity, which may be as large as a factor of ten, 
and also a model dependent factor. See App. A for this issue. There is also numerical factor 
that is derived from numerical calculations of Boltzmann equations in Ref. [39] . Here we 
collectively denote these factor as c, i.e., we write = ca^q^Jk and discuss the 

uncertainty of the maximal temperature of the Universe. Since we are mostly interested 
in the case of f j ^ <C 1, we focus on the effective termperature at the time of t ~ 
tjnax- The effective temperature T* becomes the maximal value at the time when rspiij(fc 
rrij) cii H is satisfied. Note here that there are another two model dependent factors: (i) 
total degrees of freedoms of plasma, g*, and (ii) a factor of the diffusion coefficient, = 
[See Eq. (A.5)]. As a result, the maximal temperature is given by the second line of 
Eq. (3.71) times the uncertainty factor of [(36/7r^g*)cc^/^]^/^. This dependence is rather 
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small because the time-dependence of the effective temperature has a small power of — 1/4. 
We may conclude that the uncertainty of the maximal temperature is at most a factor of 
ten. 


4 Applications 

In this section, we discuss implications of obtained results to the dynamics of scalar con¬ 
densates in the early Universe. In particular, we focus on their effects on SSB and onset of 
oscillation, which are important to determine the fate of the Universe. 

Suppose that <p has a tachyonic mass at zero-temperature and is a field responsible 
for the SSB of some symmetry. When it acquires a thermal mass larger than its zero- 
temperature tachyonic mass, it is stabilized at the origin of the potential. As the tempera¬ 
ture decreases, the thermal mass becomes comparable to its zero-temperature mass. Then 
(p starts to oscillate around the true vacuum and obtains a nonzero VEy so that the SSB 
occurs at that time. The SSB may results in formation of topological defects, such as cosmic 
strings and domain walls. These topological defects may predict some detectable signals. 
In particular, the energy density of domain walls decreases as a~^, so that they eventually 
dominate the Universe and spoil the success of the standard cosmology [19]. However, 
when the SSB occurs before inflation, these topological defects are washed out by inflation. 
Thus, it is important to determine the condition that the SSB occurs after inflation. 

Also, a scalar field can have a large expectation value during inflation due to the nega¬ 
tive Hubble induced mass term [7]. Such a scalar condensate starts to oscillate around its 
effective potential minimum when the potential force becomes comparable to the Hubble 
friction. The subsequent evolution of the Universe can crucially depend on the time when 
the scalar field starts to oscillate. For instance, if its oscillation starts earlier due to the finite 
density effects, then the scalar field tends to dissipate its energy earlier into the background 
plasma. More importantly, the baryon asymmetry of the Universe is fixed at the time when 
the AD field starts to oscillate, as we will see later. 

In the next subsection, we calculate the condition that the SSB occurs after inflation. We 
apply the result to a QCD axion model, which is well motivated in light of solution to the 
strong CP problem. Then we also apply the result to the electroweak symmetry breaking 
and comment on the DM production from the thermal plasma. Finally, we consider the 
dynamics of a flat direction in SUSY theories, especially in the context of the Affleck-Dine 
baryogenesis. 


4.1 PQ S5mametry 

In this section, we consider the following potential for the field (p: 

Vi(p) = ml\(p\^ + — (j(p\^, (4.1) 

where we include the thermal mass term (m^ ~ ol^T^ ). The SSB occurs at the effective 
temperature of f’sSB that is determined by 771^(7^^^) — 

When we require mj{T) < Av throughout the history of the Universe, we obtain an 
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(4.2) 


upper bound on the reheating temperature as 


Trh < 2 X 10'^^ GeV 


(a 


-1 


Xv 


m. 




1/2 


10^2 GeV ) Gevj 


where we assume that fj > a^, a a^, and that particles which couple to cf) (PQ-quarks) 
are not produced directly from inflaton decay. This implies that the reheating temperature 
has to be much lower than the dynamical scale v to avoid the symmetry restoration of that 
symmetry.*^^ 

Here, let us consider a QGD axion model with right-handed neutrinos [58, 59, 60, 
61] and identify cf) as the field responsible to the SSB of PQ symmetry. When the SSB 
occurs after inflation, cosmic strings form at the phase transition [62] . After the QGD phase 
transition, the non-perturbative effect associated with instantons breaks 1/(1 )pQ symmetry 
down to Z„, where n is an integer depending on models. This implies that domain walls 
form at the QGD phase transition. While these domain walls are short lived in the case 
of n = 1 [63, 64, 65, 66], they are stable and disastrous in the case of n > 2 [19, 20]. 
One of the simplest solution of this domain wall problem is that the PQ symmetry is never 
restored after inflation.This scenario requires a reheating temperature lower than the 
one derived in Eq. (4.2). QGD axion models predict a pseudo-NG boson called axion, which 
is a good candidate of DM. The abundance of axion is related to the PQ breaking scale v, so 
that the observed DM abundance determines its value [67, 68, 69]. The result is given as 


V ~8 X 10^^ GeVx n|0of, (4.3) 

where 6q is the initial phase of axion field and n is the domain wall number. This implies 
that the reheating temperature should be lower than 10^° GeV [see Eq. (4.2)]. Hence, lep- 
togenesis maybe marginally realized to explain the baryon asymmetry of the Universe [70] 
(see also Ref. [71]). 


4.2 Electroweak phase transition and DM production 


In order not to spoil the success of the BBN, reheating has to be completed before the BBN 
epoch, which means MeV. This requires that f j has to be larger than the following 

value: 


^ f y r _y" 

m3 yiMeVj yiO^^GeVj 


In this case, the effective temperature of the soft sector can be as large as 


a 


4/5^2/5 


I rnj 


2/5 


200 GeV X a^/3 


10 
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10^3 GeV j ■ 


(4.4) 


(4.5) 

(4.6) 


This is merely a necessary condition to avoid the symmetry restoration. If the field value during inflation 
is much larger than v due to the negative Hubble induced mass, the non-thermal phase transition can occur, 
which may result in formation of topological defects [54, 55, 56, 57]. To study its fate, one may also have to 
take into account of background plasma, which is involved and beyond the scope of this paper. See [37] for 
instance. 

In this case, the axion DM may predict sizable isocurvature fluctuations [21, 22, 23], which is tightly 
constrained by the observations of CMB fluctuations. This severely restricts the energy scale of inflation to 
suppress isocurvature fluctuations. 


29 










Assuming a = 10“^, we obtain T* ~ 30 GeV for the case of irij = 10^^ GeV. Even if 
-particles are generated directly from the inflaton decay, the effective temperature of the 
hard sector is at most 
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r/ ruj 


0.3 GeV X 
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10 


-27 
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(4.7) 

(4.8) 


just after inflation, where we assume t ~ 1. 

Here, let us consider thermal production of DM in a low reheating temperature. Even 
if the reheating temperature is lower than the mass of DM, it can be produced from the 
thermal plasma before the reheating is completed [28, 29] . The present energy density of 
DM from the thermal production divided by the entropy density of the Universe is calculated 
as 



(4.9) 


where the subscript “E” represents the corresponding value at the time of DM freeze-out 
(see, e.g.. Ref. [33]). To derive this, it is assumed that the maximal temperature of the 
Universe after inflation is larger than the mass of DM. However, Eqs. (4.6) and (4.8) show 
that the maximal temperature of the Universe is at most the electroweak scale for 
^(1) MeV. This means that in such a low reheating temperature the thermal production of 
DM is not efficient and Eq. (4.9) cannot be applicable to calculate the amount of DM. Even 
in this case, DM is produced from interactions between the hard and soft sector and this 
contribution is actually much more efficient to generate DM [33] . 

Next, let us identify the field cf) as the Higgs boson [14, 15, 16, 17, 18]. Then X4>- 
particles are identified as the quarks, leptons, and 51/(2)^ x U(l)y gauge bosons. In gen¬ 
eral, they are generated directly from the inflaton decay, so that the case of “direct Xij,- 
production” should be applied to this case. Eirst, the Higgs boson obtains a thermal mass 
from hard -particles just after inflation (t ~ 1) as T* with Eq. (4.8), but this maybe 
too small to restore the electroweak symmetry. Then, at the time around i = the ef¬ 
fective temperature of the soft sector reaches a maximal value and the Higgs boson obtains 
a thermal mass from -particles as T* with Eq. (4.6). This can be of order 100 GeV 
for the case of rrii = 10^^ GeV if the reheating temperature is slightly larger than 1 MeV, 
say Trh 5 MeV, because ^ y^/47t ~ 0.1 and a ~ 0.1, where and are 
the top Yukawa coupling and the strong coupling, respectively.*^^ Thus, the electroweak 
symmetry may be restored at the time around i = even if the reheating temperature is 
as low as ^(1) MeV and the inflaton mass is 10^^ GeV. When the Higgs field stays at the 
symmetric phase, the sphaleron effect, which is a 5 -I- L violating process in a finite temper¬ 
ature plasma, is turned on [72]. Our result implies that there may be an era in which the 
sphaleron effect proceeds efficiently even if the reheating temperature is ^(1) MeV. This is 
an important result when one determines the baryon asymmetry of the Universe in such a 
low reheating temperature scenario. 

Although the top quark is fermion, its effective temperature at a time around i ~ coincides with the 
one for bosons, which we calculate in this paper. 
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4.3 Affleck-Dine field 


In SUSY theories, there are many scalar fields whose potentials are absent in the limit of 
exact SUSY and within the renormalizable level. Focusing on a baryonic scalar field with 
such a fiat potential, called the AD field, baryon asymmetry can be generated by the Affleck- 
Dine mechanism [6, 7]. In this subsection, we identify (p as the AD field. 

The AD field (p has a large VEV due to a negative Hubble-induced mass during inflation 
and the infiaton oscillation dominated era. When the Hubble parameter becomes compara¬ 
ble to the curvature of the potential, it starts to oscillate around the origin of the potential. 
At the same time, it is kicked in the phase direction and generate baryon asymmetry. The 
resulting baryon-to-entropy ratio Yf, can be calculated from 

4^osc 

where we neglect 0(1) numerical factors (see, e.g., Ref. [73]). Here, 1712/2 is gravitino 
mass, is the VEV of (p at the time of beginning of its oscillation. Note that the amount 
of baryon asymmetry depends on the Hubble parameter at the beginning of oscillation 
Thus, we should include thermal effects on the potential of the AD field to determine 
as 




^3/2 

osc 


(4.10) 


H Max 

osc 


mz 


a 


4 ^ osc 


(4.11) 


where we assume \g^(p\ T [26, 27]. In the literature, they use the relation of T = 
(r^MpjH)^/"^ to calculate [see Eq. (3.59)]. However, as shown in this paper, the finite 
time scale of thermalization implies that this relation holds only after the time of t = 
so that we should check that H~^m^ > is satisfied. 

Assuming T = (Tj^MpiH)^/"^, we can rewrite Eq. (4.11) as 


Max 


m^, a 


2 


k 


2 

osc J 


We require the following condition to use the relation of T = (TpuMpiU)^^"^: 


TT < „16/5p3/5 
^osc ~ ^ ^ I 


(4.12) 


(4.13) 


Using Eq. (4.10) to eliminate \ (p\osc assuming ^ a^T^M^i/ |(/>osck we can rewrite 
the condition as 


Trh > 4 GeV X 




(4.14) 


This is easily satisfied, so that we can consistently use the relation of T = (T^MpiH)^/'^ to 
calculate the baryon asymmetry. 


Einally we comment on the case of “direct -production”. When the VEV of the AD 
field is smaller than the infiaton mass, -particles maybe produced directly from the decay 
of infiaton. In this case, the AD field may obtain a large effective mass as Eq. (4.8) just after 
inflation. Since this thermal mass easily exceeds the mass of the AD field m^, calculations 
of the baryon asymmetry may be completely changed in this case. We leave this issue for a 
future work. 
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5 Conclusion 


We have investigated the reheating process after inflation, and derived the evolutions of 
distribution function of light particles -particles), and also effective temperatures that 
describe thermal effects on scalar fields. Our results are summarized in Sec. 3.5 and are 
illustrated as the blue lines in Figs. 10 and 11. Comparing with the result calculated by 
assuming the “instantaneous thermalization”, which is widely assumed in the literature, we 
find that the maximal temperature of the Universe is overestimated in the literature, in 
particular, for the case of = FjKm ^This is mainly because emissions of 

high momentum particles are necessary for primary particles to be thermalized but such 
processes are suppressed by the LPM effect in reality. 

We have applied the result to some important phenomena: the SSB of PQ symmetry, 
the electroweak symmetry breaking, and the dynamics of an AD field. We have found that 
the maximum temperature of the Universe can be much lower than the conventional esti¬ 
mation where the “instantaneous thermalization” is assumed. In particular, we have shown 
that if the reheating temperature is as low as ^(1) MeV, then the maximum temperature 
of the Universe may be at most electroweak scale. This implies that the electroweak sym¬ 
metry may be marginally restored after inflation even for the case of such a low reheating 
temperature. The PQ symmetry might not be restored after inflation in the axion cold dark 
matter scenario, even if the reheating temperature is as large as the one required by the 
realization of thermal leptogenesis. Our results are applied to the calculation of the baryon 
asymmetry in the Affleck-Dine baryogenesis. We have justified the conventional calcula¬ 
tions performed in the literature when the VEV of the AD field is larger than the mass of 
inflaton. Also, our result implies that DM may not be able to be produced from the thermal 
plasma in a low reheating temperature scenario, contrary to the conventional studies under 
the “instantaneous thermalization” assumption. 

We have studied in detail the evolution of distribution function of radiation, and how 
particles with masses of |g^0| are produced from radiation. Thus, our results are also 
applicable to heavy particle production in the detailed process of thermalization, as partly 
done in Ref. [33]. 
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A LPM-suppressed Splitting Rate 

Here we summarize basic formulae of LPM-suppressed splitting rate in SU(N) gauge theory 
with Np-flavor fundamental fermions to understand the model dependence of splitting rate. 
We follow Ref. [51] in the following. 

Let us first start with a precise form of the collision term for the splitting process. When 
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the primaries have larger energy than the typical screening mass scale of plasma, like in 
our case, we have to take into account nearly collinear emissions. The collision term in this 
case can be expressed as 


'^split[/a] 


(2nr 

2k^ 


V. 




“ b,c J 


X (/„((c) [1 ±A(p')] [1 i/cOT] - [1 ±/„(fc)] 


<*P'<*P*(*= + P-P')rai (?'!*=•?) 


k^v 


X ifaWfbiP) [1 ±/c(pO] - [1 ±/a(fc)] [1 ±/b(p)] /c(pO) , 

(A.1) 


where represents the number of degree of freedom for species a (normalized by one real 
d.o.f.). The splitting/joining process for a be is characterized by a splitting function 
As explained intuitively in Sec. 3, we have to include destructive interference effect, 
LPM effect, in order to correctly cope with the collinear splitting process in medium. The 
required splitting function which deals with the LPM effect at the leading logarithmic order 
is given by [74, 51] 


(P; (1 - x) P, xP) = rl (P; (1 - x) p, xp) 
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(A.2) 

(A.3) 

(A.4) 


where Pj^jj^ = Min [1,(1 — x)] P, d^ is the dimension of representation for a species a, and Q 
is the quadratic Casimirs for a species a. q and g represent a fermion with the fundamental 
representation and a gauge boson of SU(A1) respectively. The diffusion coefficient is defined 
as 


Pel 




V t 
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faiP) [l±/a(P)] , 


Jp 


(A.5) 


with Vq = v^/d^ being the number of d.o.f. of a species a excluding the “color” factor 
and being the normalization of representation for a species a. The Leading Logarithmic 
enhancement factor is given by 
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where the effective temperature is defined as 


T^efi = 


T.a'^ataijaiP) [l±/a(P)] 




faip) 


(A.7) 


For instance, for a thermal-like distribution ~ Fj/p, like Eq. (3.16), the effective tempera¬ 
ture, Feff, coincides with the soft temperature, T^. 

It is instructive to reproduce the important equation [Eq. (3.15)], which is responsible 
for production of the soft population from hard primaries, and also determines 
hence To be concrete, let us assume that inflaton decays directly into SM gauge 

bosons. The Boltzmann equation responsible for production of soft particles is given by 


(271) 
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(A.9) 


Thus, the splitting rate can be parametrized as 


rspiit(k) = c'a 


9el 

~k 
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, ^/2C. 

with c = —I—LL(0;A,A,A). 




(A.10) 


The coefficient c' is given by a model dependent factor times the leading log enhancement. 
One can derive a similar formula when the inflaton decays into fermions and then it emits 
gauge bosons collinearly. Also, Eq. (3.5) can be derived from Eq. (A.l). 

In the main part of this paper, we have omitted this factor to avoid model dependent dis¬ 
cussions, but one can see that this factor results in at most ^(10) uncertainty. As discussed 
in the end of Sec. 3, the maximum temperature has a mild dependence on this uncertainty 
factor, oc Together with the numerical factor derived from numerical calculations of 
Boltzmann equation, we conservatively conclude that the uncertainty of maximum temper¬ 
ature is at most a factor of ten. 


B Effective Temperature for 2 '(/,-particles 

In contrast to -particles, -particles can have sizable masses of \g^4>\ due to the large 
expectation value of the (/)-condensate. The distribution in momentum space can be differ¬ 
ent from that of -particles and hence separate discussions are required. In this Appendix, 
we show the results of T* . 


B.l <t <t 


-mi 


soft 


If the effective mass of -particles are smaller than the screening mass, then the effective 
temperature for -particles coincide with that for -particles: 
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for \g^(p\ < m^. Here the subscript “indir” indicates -particles which are not produced 
directly from inflaton decay. For the case of < |g^0|, the effective temperature squared 
is evaluated as*^® 


= Max 
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By using these equations, one can show that, for a ^fcj^ax ~ \g 4 ><P \ with < a°, the contri¬ 
bution of abundant -particles is always subdominant compared with that of the running 


coupling constant. Even in the case of Fj > a°, its effect is roughly the same order with 
that of the running coupling constant for |g^(/)| > Hence, one can omit the decay 

of massive -particles in order to show that the effect on the effective potential of (p from 
-particles is subdominant for the large field value regime. 


However, if -particles are produced directly from inflaton decay, it has a sizable con¬ 
tribution to the effective temperature, which tend to dominate for a large field value [see 
Eq. (3.40)]. Therefore, we should take into account the decay of -particles to obtain 
more realistic predictions at least in that case: 



dir 


= arV^T^ 




(B.6) 


where the subscript “dir” indicates contribution from -particles which are directly pro¬ 
duced from inflaton decay. Hence, one should compare it with the log potential from the 
running coupling constant. The effective mass squared for the large field value regime, 
|g^(/)| > should be given by 
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If the effective mass of -particles are smaller than the screening mass, then the effective 
temperature for -particles coincide with that for -particles: 



indir 


(B.8) 


Here we neglect t-channel contributions in the soft sector, because they are always subdominant to 
calculate the thermal potential. 


35 




































for \g^(p\ < m^. Here the subscript “indir” indicates -particles which are not produced 
directly from inflaton decay. For the case of rUj < 01, we obtain the effective temperature 

squared as 
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where represents the interval boundary the hard and soft sectors. Here we include the 
contribution coming from t-channel scatterings between the hard and soft sector. Again, 
one can show that the effective potential from -particles are always subdominant for the 
case of I01 > with < a~^. And also even in the case of > a~^, it is roughly 

the same order of that from the running coupling constant. Thus, we can omit the decay 
of -particles in order to show that it is always subdominant compared with that from the 
running coupling constant for the large field regime. 

However, if the inflaton directly decay into -particles, then the effective temperature 
from direct inflaton decay tends to dominate for the large field value, |g^0| > 

Hence, one should compare it with that from the running coupling constant: 


Max 
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(B.13) 


where |(jir is given by Eq. (B.6) [See also Eq. (3.54)]. 

For the small field value regime of |g<^^| < T^,, one can show that the soft -particles always dominate 
over that of hard ones produced via direct inflaton decay. 
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^•3 tmax < t < tRH 


obtain the effective temperature squared 


for \g^4)\< r* 

for r* < |g^(/)| < ( 3 ^ 4 ) 

where the subscript “indir” indicates -particles which are not produced directly from 
inflaton decay. At that time, one can show that the contribution from abundant -particles 
is always subdominant compared with that from the running coupling constant for the large 
field value regime |g^(/)| > a~^T^. Thus, we can omit the decay of -particles in order to 
show that it is always subdominant compared with that from the running coupling constant 
for the large field regime. 

However, if the inflaton directly decay into -particles, then its contribution tends to 
dominate for the large field value regime, |g^0| > Hence, one should compare it 

with that from the running coupling constant: 


Calculating the distribution of -particles. 



Max 


T4 




02’ *.X^ 

dir 


(B.15) 


where is given by [See also Eq. (3.61)] 




Trt ^m?Min 


dir 


_^split(^/)^ ^decay^ 


decay_ 


(B.16) 


Here note that the hard particles with the momentum of soon breaks up within the time 
scale of Tgpiit, and hence l/[rspijt(m|)t] should be multiplied. 
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